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Locomotion of bacteria by actin polymerization, and in vitro motion of spherical beads coated 
with a protein catalyzing polymerization, are examples of active motility. Starting from a simple 
model of forces locally normal to the surface of a bead, we construct a phenomenological equation 
for its motion. The singularities at a continuous transition between moving and stationary beads 
are shown to be related to the symmetries of its shape. Universal features of the phase behavior 
are calculated analytically and confirmed by simulations. Fluctuations in velocity are shown to be 
generically non-Maxwellian and correlated to the shape of the bead. 
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Motility, the capacity of organisms for independent 
motion, has long fascinated physicists as well as biol- 
ogists. For example, the reversible character of low 
Reynolds number flows [1] requires ingenious methods for 
swimming motion [2]. Another strategy for movement, 
the subject of this paper, is propulsion via forces gener- 
ated by polymerization of actin on an organism's surface. 
Polymerizing actin filaments are involved in, among other 
things, the locomotion of certain bacteria, such as the 
pathogenic Listeria monocytogenes. Upon entering the 
cytoplasm of a host, a Listeria "hijacks" the cell's actin 
machinery; first a roughly symmetrical cloud of polymer- 
izing actin filaments envelops the bacterium, and then an 
asymmetric "comet tail" of actin spontaneously develops 
that propels the bacterium forward [3] . Motility by actin 
polymerization has been the subject of several analytical 
and model studies [4-8]. 

A number of in vitro experiments have been performed 
on artificial systems designed to mimic Listeria motil- 
ity under simplified and more easily controlled circum- 
stances [9-13]. In these experiments, the bacteria are 
replaced with spherical beads coated with a protein that 
catalyzes actin polymerization. When these functional- 
ized beads are placed in cell extracts or a mixture of 
purified proteins, phenomenology similar to that of Lis- 
teria motility is observed. A spherically symmetric cloud 
of cross-linked actin filaments forms around each bead. 
Then, depending on the various parameters of the exper- 
iment, some beads are observed to spontaneously break 
free of their actin clouds and move away, propelled by an 
asymmetric actin comet tail. 

Several theories have been put forward to explain the 
mechanism behind this spontaneous motility, van Oude- 
naarden and Theriot model the actin filaments surround- 
ing the bead as elastic Brownian ratchets [14]. The 
stochastic growth of these filaments is coupled due to the 
presence of the bead; when one filament grows, it gives 
the bead a small push, making room for neighboring fil- 
aments to grow. This cooperativity generates positive 
feedback and the possibility of an avalanching symmetry- 
breaking event. In another model, Noireaux et al have 
proposed that the actin cloud surrounding the bead be 



treated as an elastic gel [15]. Polymerization occurs at 
the bead surface, which expands the gel and stretches it. 
At a certain point, stress at the outer surface ruptures 
the gel, and the bead is free to escape. Mogilner and 
Oster have suggested a third alternative in which move- 
ment occurs because of the failure of cross-links within 
the actin network [8] . 

Our aim is not to reconcile these existing theories, or 
to provide a comprehensive description of the observed 
phenomena. Instead, we investigate the behavior of an 
activated rigid bead within a framework that is sufficiently 
general to include many possible microscopic models. In 
the spirit of the Landau theory of magnetism, we treat 
the actin filaments surrounding the bead as a continu- 
ous effective field that generates a normal force density 
at each point on the bead surface. The evolution of this 
order parameter field is coupled to the velocity of the 
bead, allowing the possibility of positive feedback as in 
the above Brownian ratchet model [14]. Our approach 
neglects many aspects of the system, but can provide an- 
swers regarding generic features: What are the possible 
singularities that can accompany the transition between 
stationary and moving states, and how do they depend 
on the shape of the bead? Regarding the active bead as a 
stochastic system, are there characteristics of its fluctu- 
ations that set it apart from a passive (Brownian) bead? 

We consider a bead of fixed shape, with unit normal 
vector n(r) at a point r on its surface. The polymer- 
ization of actin (or any other actively energy-consuming 
process for that matter) is assumed to exert a force, lo- 
cally directed along the normal, at each point on the bead 
surface. The force per unit area at time t is indicated by 
a scalar field g(r, t), which could for example be a func- 
tion of the actin density at that point. The net force due 
to such activity, acting on the center of mass of the bead, 
is F a (t) = J dS g(r,t) n(r), the integration being over 
the entire bead surface. In response to force, the bead 
moves at an instantaneous velocity v(t). In a viscous 
medium the components of the force and velocity are lin- 
early related by v a = ^. a pFp with a, (5 = 1,2,3, where 
the mobility tensor fi a p depends on the shape of the ob- 
ject. For functionalized beads propelled by a comet tail, 
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drag from tethered filaments in fact far exceeds viscous 
drag [11]. Whether this would be true for a stationary 
bead in an actin cloud is unknown. However, we are rel- 
atively unconcerned about the exact nature of the drag 
force: For small enough velocities, it should be possible 
to make a linear expansion in force, with higher-order 
corrections appearing at higher speeds. 

To complete the description of the dynamics, we need 
to specify the time evolution of g. We assume that the 
rate of change of g(r, t) is a local function of <?, and that it 
depends on the net velocity v. The dependence on v cou- 
ples the field g at different positions on the bead and also 
provides the possibility of positive feedback (e.g. when 
polymerization is decreased on the front, and enhanced 
on the rear of a moving bead). The growth rate is then 
expanded in powers of g and v, resulting in 



dg(r,t) 
dt 



= $ [<?(r), v] = — g + av • n — g 2 + bgv ■ n 
+dv • v + e (v • n) 2 — eg 3 + ■ ■ ■ + T](r, t), (1) 



where we have included all terms to second order, and the 
simplest cubic term. Note that the vector v contributes 



through scalar forms i>j 
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an expansion is presumably valid close to a continuous 
transition in which both g and v are small, and needs 
to be sclf-consistcntly justified. We have rescaled g and 
t so that the coefficients of g and g 2 are both — f . The 
sign of the linear term is negative so that a uniform actin 
distribution (set to g = without loss of generality) is 
stable in the absence of any velocity coupling, while the 
sign of the g 2 term is unimportant since we may always 
consider the dynamics of — g rather than g. In addition 
to rescaling g and t, we may also rescale lengths; a con- 
venient choice is to set the surface area of the bead to 
unity. Through appropriate redefinitions, our previous 
formulas for the force F a (t) and bead velocity v remain 
unchanged. For analysis and numerical simulations of 
Eq. (1) we typically set d = e = ■ ■ ■ = 0, and c > for 
stability. The final term in Eq. (1) allows for a stochastic 
noise rj(r, t). 

In the absence of an external force, Eq. (1), along with 
the formula for velocity 



v a (t) = n a pFp{t) = a af) J dS g(r,t) n p (r), 



(2) 



fully specifies our model. We first analyze the case of a 
spherical bead, and then generalize its lessons to more 
complicated shapes. Due to its symmetry, the velocity 
and force are parallel for a sphere, and fi^p = n8 a p- 
In the absence of noise 77, there is a trivial solution of 
g(r, t) = v(t) = 0, which corresponds to a static, uniform 
actin distribution and a motionless bead. This solution 
is linearly stable for a < 3/fi, but unstable to dipolar 
fluctuations for a > 3//i. In the latter case, the unsta- 
ble actin fluctuations grow and saturate. Thus Eq. (1) 



allows for a symmetry-breaking transition from a bead 
at rest (v = 0) to a bead in motion (v = const) as the 
parameter a is varied past a critical value a c = 3/fj,. This 
parameter controls the strength of the positive feedback 
in the coupling of the bead velocity to the actin field. 

Near this transition we expand g in spherical harmon- 
ics with g(r, t) — » g(Q, t) = J2i m 9tm{t)Y im (Q), where 
represents the solid angle coordinates 9 and <f>. In this 
basis Eq. (1) becomes 



dge 



dt 



(-1 + ^~<^i) 9tm + nonlinear terms, (3) 



where the nonlinear terms couple different gem's 
with coefficients of the form J dSlYe imi Ye 2m2 Yg m , 
J dftYe imi Ye 2m2 Ye 3m3 Y e * m , From the linear terms 

in Eq. (3) we see that, near the transition with e = 
(i(a — a c )/3, |e| <C 1 and for small actin fluctuations, 
the dipolar (I = 1) modes evolve on a slow time scale 



while the other modes evolve on a fast time 



scale t cx 0(1). Due to their instability (or near instabil- 
ity) we also expect the dipolar modes to be larger than 
the others. We therefore treat the fast t 7^ 1 modes as 
adiabatically slaved to the slow I = 1 modes and self- 
consistently solve for their amplitudes as functions of the 
<7i m 's. Substitution then yields effective equations of mo- 
tion for the slow modes, and by a simple transformation, 
the bead velocity [20], as 
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with v 3 = w 2 v, etc. Thus, near the transition the bead 
velocity v obeys a Ginzburg-Landau equation consistent 
with spherical symmetry. When the coefficient of the cu- 
bic term is negative, the transition to nonzero v is con- 
tinuous with the steady state velocity v* going to zero as 
e — > + . In the vicinity of this continuous transition in 
parameter space, our analysis above is self-consistently 
justified and Eq. (4) should provide an accurate descrip- 
tion of the dynamics. The fifth-order coefficient u, al- 
though not explicitly presented above, is negative in this 
region, providing overall stability. 

The phase diagram and critical behavior associated 
with the steady-state solutions v* of Eq. (4) are fa- 
miliar from mean-field theory in equilibrium statisti- 
cal mechanics. The cubic coefficient is negative when 
feci < b < fe C 2, where b c \ and fe C 2 are the two roots 
of (fib/'i - l)Gufe/3 - 2) = c (recall that c > 0). In 
this regime, the transition is continuous with v* cx e 1 / 2 . 
When the cubic coefficient is positive, the transition is 
discontinuous. Tricritical points, at which v* cx e 1 / 4 , oc- 
cur when 6 = b c \ or fe = fe C 2 precisely. Along the line 
a = a c or e = 0, v* cx |fe — fecil 1 ^ 2 for i = 1, 2 as b varies 
outside of the interval [fe c i , fe c2 ] . Finally, in the vicinity 
of the first-order transition with 6 > fe C 2 or 6 < fe c i, the 
v = and v — const states coexist, both being locally 
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FIG. 1: Schematic phase diagram for a spherical bead in the 
vicinity of continuous transitions for a — a c and b c \ < b < 
b C 2- In the grey region both stationary (v = 0) and moving 
(v = const) states are locally stable. 

stable. In the a-direction, this region extends over the 
range a* < a < a c where the "limit of metastability" a* 
has the b dependence a c — a* oc (b — b c ) 2 for b < b c \ or 
b > b C 2- The schematic phase diagram for the motion of 
the bead is shown in Fig. 1. 

To confirm the analytical results, we numerically sim- 
ulated Eq. (1) with a fourth-order Runge-Kutta scheme. 
For integrations over the surface of the spherical bead, 
we adapted routines from the NAG library [16]. We con- 
firmed the phase diagram in Fig. 1, as well as the various 
critical exponents. As an example of the latter, Fig. 2 
illustrates the time evolution of velocity for a > a c = 1. 
Starting from an initial (unstable) stationary state, any 
small fluctuation takes the sphere into a moving phase. 
The characteristic time scale for this change of state di- 
verges as e _1 cx l/(ct — a c ), while the saturation velocity 
scales as e 1 / 2 oc y/a — a c . As indicated in Fig. 2, the ve- 
locity evolution curves obtained for a number of different 
a can be collapsed together using the above scalings of 
velocity and time. 

Let us now consider a rigid bead of arbitrary shape. 
The acceleration of the bead can be obtained by taking a 
time derivative of Eq. (2) and substituting from Eq. (1). 
The analog of Eq. (4) for the sphere is then obtained as 

dv 

—jjT = (S a p + a [i ai njnp ) vp H h /«(*)• (5) 

At the linear order, properties of the shape are encoded 
in the tensor n a np = J dS n a (r)np(r). The stochastic 
force f a (t) partly originates in the noise in Eq. (1), and 
correlations amongst its components also depend on the 
shape through the matrices /i a/ g and n a np. The linear 
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FIG. 2: Rescaled bead speed versus rescaled time for /i = 3, 
c = 0.2, b — 1.5 (bd ~ 0.83 and b C 2 ~ 2.17), and various 
values of a > a c = 1. The bead begins (very nearly) at rest 
at t = with a small random actin distribution. The curves 
are nearly perfectly superimposed. 



stability of a stationary bead is determined by the eigen- 
values of the matrix fi ■ nfi. When a\i > 1, the 
bead is unstable to motion in the corresponding eigendi- 
rection. The largest eigenvalue indicates the direction 
that initially becomes unstable. Of course, as in Fig. 1, 
it is possible to have coexistence of moving and station- 
ary states before the onset of this instability. In fact, the 
analysis so far does not rule out coexistence of several 
states moving along different directions [21]. 

It is interesting to note that eigenvalues of rTn are 
larger along directions where the shape displays a big- 
ger cross-section. Thus, ignoring variations in mobility, 
a pancake shape prefers to move along its axis (breaking 
a two- fold symmetry), while a cigar shape moves per- 
pendicular to the axis (breaking a degeneracy in angle). 
This is opposite to what happens for passive ellipsoids 
in a fluid where hydrodynamic effects favor motion along 
the slender axis. 

To see what happens when the stationary bead is 
linearly unstable, we need to examine the higher-order 
terms in Eq. (5) [22]. // the transition is continuous, 
its singularities are determined by the next higher-order 
term in the expansion. For the sphere, symmetry con- 
siderations rule out a quadratic term; the cubic term 
leads to the singularities discussed earlier. However, 
quadratic terms vpv~f may appear in Eq. (5) with coeffi- 
cients C a f3j that are related to the shape by n„n ( 3?J 7 = 
J dS n a (r)n/j(r)n 7 (r). If such terms are present, the 
velocity will vanish on approaching the transition as e, 
rather than y^e. To distinguish the two universality 
classes along a particular direction, we merely need to 
ask if the opposite direction is equivalent. Thus the ve- 
locity of a cigar-shape will vanish as y/e, while that of 
an arrowhead goes to zero linearly. In the language of 
dynamical systems, the transition at e = is a pitchfork 
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bifurcation in the former case and a transcritical bifur- 
cation in the latter [17]. 

Let us now consider the effect of the noise f a (t) in 
Eq. (5). At a coarse level this resembles a Langevin 
equation for the velocity v. However, since this dynamic 
equation is not subject to fluctuation-dissipation con- 
straints, the resulting probability distributions need not 
be those familiar from equilibrium statistical physics. For 
example, let us consider an ellipsoidal shape in the regime 
where stationary motion is stable, and higher-order terms 
in Eq. (5) can be ignored. For a passive (Brownian) ellip- 
soid, the velocity fluctuations follow a Maxwell probabil- 
ity distribution. In this case the center of mass velocity is 
equally likely to point in any direction, irrespective of the 
orientation of the ellipsoid. For active fluctuations (e.g. 
due to actin polymerization) of the same bead, the prob- 
ability distribution of velocities can be extracted from 
Eq. (5). Linear truncation leads to a Gaussian probabil- 
ity distribution, P(v) cx cxp (— v a M a pVfj), but there is 
no a priori reason for the matrix M a p to be proportional 
to S a fj. Thus generically the velocity fluctuations of the 
active bead are non-Maxwellian, and correlated with the 
shape of the bead. 

In summary, motivated by experiments on actin-based 
motility, we have characterized some of the phenomena 
that distinguish such an active system from a passive 
Brownian particle. We start with a highly simplified mi- 
croscopic model whose central ingredients are a force that 
is (a) locally normal to the surface of a rigid bead, and 
(b ) influenced by the velocity of the bead in the surround- 
ing medium. These assumptions lead to a macroscopic 
equation for the bead velocity whose coefficients reflect 
the symmetries in its shape. (The reliance on symmetry 
should increase the validity of the form of this equation 
beyond our particular microscopic model.) We note that 
Eq. (1) for the dynamics of the actin distribution is sim- 
ilar to equations that model pattern formation [18] or 
"self-organization" [19] in various contexts. As in our 
case, analysis of these systems near threshold can often 
be performed by eliminating fast modes. 

With the macroscopic equation for the bead velocity 
we can analyze the singularities that accompany a con- 
tinuous transition between stationary and moving states. 
Depending on symmetries along the direction of move- 
ment, the velocity vanishes with a square-root or linear 
singularity. We also examined the fluctuations in velocity 
of an active bead. We find that, generically, the proba- 
bility distribution is not Maxwellian in form, and is de- 
pendent upon the orientation of the bead. It would be 
interesting to experimentally determine this distribution 
though single-molecule experiments on active beads. In 
this work we focused only on the translational motion of 
a freely moving bead. Rotations of a bead, as well as 
its motion in a trap, are likely to bring in additional de- 
grees of freedom (position and angular coordinates) that 
are coupled to the velocity of the bead. In future work, 



we would like to incorporate such effects to form a more 
complete picture of the repertoire of motions available to 
active beads. 
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